(N 



o 

(N 



c^ 



NONPERTURBATIVE EVALUATION OF THE DIFFUSION 
RATE IN FIELD THEORY AT HIGH TEMPERATURES 

Alexander Bochkarev "[[] and Philippe de Forcrand^ 

"'Institute for Theoretical Physics, University of Bern 
Sidlerstrasse 5, CH-3012 Bern, Switzerland. 

HPS, ETH Zurich, CH-8092 Zurich, Switzerland 



Abstract 



> 

(N 

^^ \ Kramer's approach to the rate of the thermally activated escape from 

a metastable state is extended to field theory. Diffusion rate in the 1+1- 
(N ■ 

^^ \ dimensional Sine-Gordon model as a function of temperature and friction coef- 

-)— » . 

^ ' ficient is evaluated numerically by solving the Langevin equation in real time. 

Oh, a clear crossover from the semiclassical to the high temperature domain is ob- 

i^ ' served. The temperature behaviour of the diffusion rate allows one to determine 

- ^-, . the kink mass which is found equal to the corresponding classical value. The 

rS 
Vh ' Kramer's predictions for the dependence on viscosity are qualitatively valid in 



this multidimensional case. In the limit of vanishing friction the diffusion rate 
is shown to coincide with the one obtained from the direct measurements of 
the conventional classical real-time Green function at finite temperature. 
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1 Introduction 

Statistical systems with finite energy barriers separating different domains of the 
phase space may exhibit metastable states. Evaluation of the rate of escape from 
those states is an important dynamical problem. The simple example of a statistical 
system where one degree of freedom is exposed to the potential illustrated by Fig.l 
was first considered by Kramers |jl]]. At temperatures : 

Wo < T < E5 (1.1) 

where ujq is the scale of quantum fiuctuations and Es is the barrier height, the classical 
thermodynamical fiuctuations contribute to the escape rate F ~ exp{—Es/T) of the 
system initially localized in the well. Kramers |I[ suggested to determine the escape 
rate by means of the classical Langevin equation: 

dU 
(f) = p/M, p = -— - ^p + 7]{t) (1.2) 

{V{t)vm = 2T^M5{t) (1.3) 

where 7 is a friction coefficient, introduced as an input parameter in this phenomeno- 
logical description and M corresponds to the mass of the would be Brownian parti- 
cle. 7 determines how strong is the coupling of the system to the heat-bath, repre- 
sented by the white noise r] , and therefore how fast the system reaches equilibrium. 
The normalization of the random force (|1.3| ) (the Einstein relation) comes from the 
fiuctuation-dissipation theorem. In the case of one degree of freedom the escape 
rate r(T, 7) has been obtained analytically in the domain of very small 7 and 
for moderate-to-large friction 7 ^ u;_ , where u^ is the negative eigenmode of the 
fiuctuations near the stationary point (p = (pg. In quantum field theory one often 
deals with an analogous problem of penetration through some energy barrier. One 
very important example is anomalous fermion (or axial charge) number violation 
in gauge theories with nontrivial structure of the ground state 0. Here the static 
energy barrier separates different classical vacua with definite integer Chern-Simons 



numbers [Q. Penetration through this energy barrier leads to the dissipation of the 
fermionic number. The system at high temperature is expected to exhibit thermal 
activation behaviour of the rate of anomalous fermion number non-conservation [^. 
The problem specific for the field theory is that the energy barrier - analogue of U 
in Eq. (|1.2|) - is not known explicitly. This energy barrier is in the multidimensional 
configurational space. Only the static field configuration, corresponding to the top of 
that barrier is known from topological considerations [^] or explicitly in some models 

E- 

In previous works |Q , we suggested to follow Kramers in studying the energy 
barrier in field theory. In [^ we introduced microscopical classical Langevin equation 
in real-time in the 1-1-1 dimensional Abelian-Higgs theory and observed the Brownian 
motion of the Chern-Simons variable between topologically distinct classical vacua. 
Our measurements of the diffusion rate revealed its thermal activation dependence 
on the temperature in the domain (|1 . 1| ) . Although this approach is applicable in the 
interesting domain of high temperatures, where the semiclassical approximation is 
not valid, it was difficult to get data there, because it required rather big lattices. 
Meanwhile, it was argued in |]TU[ that our results were difficult to interpret, since they 



were based on the first-order Langevin equation, derived normally for large friction 
coefficients, while, in fact, 7 = 1 had been used. 

In the present paper we deal with a simpler model, the Sine-Gordon field theory 
in 1 -|- 1 dimensions, which allows one to obtain not only accurate quantitative results 
but also to check qualitative conclusions valid beyond this particular model. In fact, 
this model is very similar to the Abelian Higgs theory. On the classical level the Sine- 
Gordon theory has an infinite number of degenerate vacua with finite energy barrier 
separating them. At nonzero temperatures one observes a random walk between those 
vacua. At temperatures ( p.. ID this random process goes through the formation of kink- 
antikink pairs, so the rate of the process (diffusion rate) is sensitive to the density 
of kinks. Previous studies of the density of kink gas as a function of temperature 
in the 1-1-1 dimensional scalar field theory with the double well potential revealed 



an intriguing fact : the thermal activation behaviour of the density of kinks n ~ 
exp{—E^-^/T) involved an effective kink mass -E^ "which was found to be 20 -i- 30% 
smaller than the classical value E"^ p| , |jTO| . One source of uncertainties in measuring 
the density of kinks was the adhoc criterion used in counting the number of kinks in a 
given field configuration. This criterion was suggested in [|l3l. One advantage of the 
Sine-Gordon model is the possibility to avoid this criterion problem. The quantity 
we measure is very well defined: it is the diffusion parameter of the average field. 
We measure the diffusion rate in the semiclassical domain ( |1 . 1| ) and beyond. The 
measurements allow us to extract the value of the kink mass E'j^, which we find to 
coincide with the classical kink mass E'j^. We therefore conclude that the previously 
observed discrepancy was an artifact of the criterion used. 

We also measure the viscosity dependence of the rate by means of the second order 
Langevin equation. We find a remarkable coincidence of this dependence with the 
one predicted by the first order Langevin equation down to very small values of the 
friction coefficient. This implies that one can use the first order Langevin approach 
in field theory, in particular, for 7 ~ 0(1). 

In the limit of vanishing friction we find a nonzero diffusion rate. Since in this 



limit the Langevin equations (|L2| ) become the Hamiltonian ones, we perform a direct 
numerical measurement of the classical real-time two-point Green function, which 
describes the diffusion, as a Gibbs's average. These ensemble averages give the same 
thermal activation behaviour of the rate as obtained in the Langevin measurements 
and, what is less trivial, they prove to be equal in magnitude to the Langevin mea- 
surements in the limit of vanishing friction. This establishes an important equivalence 
between the Gibbs ensemble measurements of the real time classical Green functions 
at finite temperature and the corresponding Langevin measurements done in the 
asymptotic domain of very small friction. 

In the next section we explain the relevance of simulations of the classical field 
theories to the behaviour of quantum field theories at high temperatures. We empha- 
size a particular relation between the lattice spacing of the classical systems under 



consideration and the temperature. In section 3 we briefly introduce the 1 + 1 di- 
mensional Sine-Gordon model at finite temperature. Then we describe the first- 
and second-order Langevin equations. Sections 4 and 5 contain the results of our 
numerical simulations. 

2 Classical systems - for quantum theories 

It is well known that the partition function oi a. D + 1-dimensional quantum field the- 
ory in the limit of high temperature may be obtained by means of the corresponding 
D-dimensional classical field theory. To see how this emerges consider, for instance, 
a scalar field theory in D + 1 - dimensions, defined by the action of the general form: 

S = jd^'+'x []^{dct>f + ^02 + Af/.J (2.1) 

The partition function of the corresponding statistical system is given by functional 
integral with the following Matsubara action: 



Sm = I dx,J d^x Uij-d,<l>)' + 2(V0)' + ^0' + XU,rA 



(2.2) 



where (3 = 1/T is inverse temperature, h is Planck constant and periodic boundary 
conditions are imposed on the field in the imaginary time direction. Eq. (|2.2|) implies 
that in any of the limits : 

n^O , or T ^oo (2.3) 

the static X4-independent field configurations dominate in the functional integral. If 
one ignores for a moment ultraviolet divergences and performs the limit ( |2.3|) formally, 
the contribution of the static modes is determined by the action 

S^ff = pjd''x f^(V0)2 + ^4>' + Af/,J (2.4) 

This action does not mention Planck constant. The functional integral with the 
action ( |2.4| ) determines the partition function of the corresponding classical system. 
To be more accurate, one should keep in mind that nonstatic modes do not decouple 



in the ultraviolet divergent diagrams, which means that the action ( |2.4| ) involves 
renormalized Largangian parameters - running coupling constants, normalized on 
the temperature. The renormalization effects are the memory about the quantum 
nature of the original theory ( p.2| ). To obtain the effective action ( p.4|) is, in general, 
a separate problem fll]]. In some simple cases in low dimensions the action (|2.4| ) is 
simply a static Hamiltonian. 

Let us see explicitly how the reduction takes place, for example, in the case of weak 
interaction (A <^ 1). The free energy of a gas of particles of mass m, corresponding 
to (p.l|) , in the one- loop approximation reads as: 

dk^ 1 



Tvj:f 



(27r)^ 2 



-log^^ + ^,^ (2.5) 



where u;„ = 27TnT,n = 0, ±1, ±2, ... are Matsubara frequencies, Uk = \rm? + W and 
V is D-dimensional volume. The temperature dependent part 

F = TVj j^log{l-exp{-f3u,)) (2.6) 

becomes in the limit T ^ m : 

r dk^ 

F = -TV J j^ (3knsim (2.7) 

where ub is Bose distribution so that 

F = -kVT^+^ (2.8) 

with 

«^ = ^r(D + i)c(z^ + i) (2.9) 



where Q = 2'k^/'^/T{D/2) is a surface of the D - dimensional sphere |T^ . 

If we want to obtain the same expression by means of the effective theory (|2.4|) we 

have to perform the limit /? — ;> in the integrand of Eq. ( |2.7|) , which gives divergent 

factor: 

r dk^ 
^--TVJ^y, (2.10) 



This divergent integral counts the number of degrees of freedom of the classical field 
theory (the Ray leigh- Jeans divergency). Thus the partition function of the classical 
field theory is something which needs to be defined. Regularization is necessary for 
this end. 

The vacuum energy of the effective theory (|2.4| ). 



TVl^^logGj, (2.11) 



(27r)^ 2 
is determined by the regularised propagator 

GA = l/{ul + u;t/A') (2.12) 

For the temperature dependent cut-off : 

A = cT, cocO(l) (2.13) 



one recovers the result Eq.( ^.8| ) with c determined by k. 

This simple exercise illustrates one general and important statement. The effective 
D - dimensional classical field theory ( |2.4| ), which serves for the calculation of the high 
temperature limit of the corresponding D+1- dimensional quantum field theory, has a 
physical cut-off, given by Eq. ( p.l3| ) . If we assume lattice regularization for the original 
theory (|2.2|) , we conclude that the continuum limit is performed simultaneously with 
the limit (|2.3|) . One should perform thermodjTiamical limit in the functional integral 
with the effective action (|2.4| ), but not the continuum one. It is essential that the 
lattice spacing a"^ oc A here is the physical cut-off, unambiguously fixed by the 
temperature: 

a-^(xT/n (2.14) 

In this sense the original quantum field theory may be studied at high temperatures 
by simulating the classical functional integral with the regularized action (|2.4|) and 
a temperature dependent cutoff (|2.14| ) |[I3[| . As soon as nontrivial lattice spacing 
dependence is observed in a study of the classical systems at high temperatures the 
relation ( p. 14 ) is to be taken into account. 
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3 Kramers approach in field theory 

We consider the Sine-Gordon model in 1 + 1 dimensions defined by the action: 

S = J d^x (^{d4>y - "^ {1 - cos[^/X4>{x)])\ (3.1) 

On the classical level the theory has an infinite number of degenerate vacua with 
(j) = 27Tn, n = 0, ±1, ±2, ... In the quantum theory the mass m determines the scale 
of quantum fiuctuations and the self-coupling constant A is bounded [jl4|: A < Sir. 
The well-known static solution to the classical equations of motion interpolating 
between different classical vacua <Pk = ^ arctan (e*""^) has energy 

Ek = 8m/A (3.2) 

which is parametrically large, in the sense that A is a parameter of the semiclassical 
(loop) expansion. If m is chosen as a unit of measure, the Lagrangian in Eq. (|3.1| ) 
may be rescaled to the one which has no free parameters: 

^ = \ h'"' ^2^^"^^' ~^^~ cos[0(x)])) (3.3) 

The Sine-Gordon model provides one with a good opportunity to study thermal 
activation phenomena in field theory. The potential energy is completely bounded at 
infinity for the zero mode (f> = jjQdxip{x), which allows for the Brownian motion 
of (f at finite temperature. Any initial state localized around one classical vacuum is 
far from equilibrium. In the case of periodic boundary conditions fiuctuations of the 
field between different vacua proceed via formation of a pair of kink and antikink and 
their separation from each other. Creation and annihilation of kink-antikink pairs 
is essentially a random walk between different classical vacua. Thus the diffusion 
rate of if is sensitive to the kink mass, which determines the height of the static 
energy barrier between different vacua. For periodic boundary conditions kinks and 
antikinks appear in pairs, so one expects in the semiclassical domain of temperatures 

TpBc ~ exp{-2Ek/T) (3.4) 



In the case of free boundary conditions a configuration with a single kink may 
appear as a result of the thermal fluctuations. The kink's motion through space 
brings the system from the neighborhood of one classical vacuum to the vicinity of 
another vacuum. Therefore for free boundary conditions we expect: 

TpBC ~ exp{-Ek/T) (3.5) 

Following Kramers, we describe evolution of some initial non-equilibrium state 
localized around one classical vacuum by solving the real time Langevin equations. 

Let us discretize the system of size L = Na , where a is lattice spacing. The 
Hamiltonian, corresponding to Eq.( p.3| ) is : 

H = K + U (3.6) 

^ = E ^ 4 iPnf (3.7) 



n=i 2 a^ 



^ = E — I 2 ('^"+1 - '^")' /«' + !- cos[(t>n] I (3.8) 

where D is dimensionality of space, D = Im our case. To generate a system with the 
density matrix p = exp (—H/Tj , one naturally employs the second order Langevin 
equations of the form: 



a 



-^Pnit) (3.9) 



Pn = ^ { {<f>n-l + </>n+l " 2(/)„) /tt^ - sin((/)„)| - 7p„ + ?]„ (3.10) 

which is a straightforward generalization of Eqs.( |L2|) , (|1.3| ) to the case of many 
degrees of freedom. Eqs.( |3.6| )-( |3.1U| ) indicate that formally we deal with a classical 



gas of particles of mass: 

M = a^/X (3.11) 

from which one concludes that the normalization of the random force, fixed by the 
Einstein relation (see Eq.( |1.3|) ), is : (^n(^)^m(O)) = 2T'ya^6(t)6nm/X . The average 
kinetic energy is independent of the effective mass M according to the usual formulae 
of quantum statistics < K > = NT/2. 
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The form of the rescaled action Eq. ( p.3| ) iinphes that the classical dynamics does 
not depend on A. To see it explicitly rescale the momentum: p = Xp. Then the 
density matrix looks like : 

p = e^v{-Heff/eeff) (3.12) 

Oeff = ^ (3.13) 

^ r 1 1 1 

Heff = E I ^(Pnf + - (0n+l - <Pnf / a^ + 1 " COs[<Prr] j (3.14) 

The corresponding Langevin equations indeed do not contain the parameter of the 
semiclassical expansion A explicitly: 

(0„(2 + l) - (Pn{^))/e = Pni^ (3.15) 

(Pnii + 1) - Pnii)) /£ = (0n.-l + <Pn+l " 20„) / O^ 

- sin(0„) - 7p„ + r]n{i) (3.16) 

where : 

{Vnii)VmiJ)) = —SnmSij (3.17) 

and we have discretized Langevin time by the amount of e. Physical quantities are 
obtained in the limit e — >■ 0. Notice that the normalization of the white noise ( p. 171) 



is different from that in |T0| 



Eqns.( prT^ )-( ^.14| ) define the effective classical system we are going to simulate. It 
incorporates the lattice spacing a as an input parameter, determined in the underlying 
quantum theory by Eq.( 2.14|) . The kinetic energy of the effective system ( 3.12 )-( p.l4|) 



averaged over the evolution ( |3.16| ) may be evaluated to check the effective temperature 

{K^jf) = ^NO^ff (3.18) 

One may, of course, consider the first-order Langevin equation instead of the 
second-order one. In the Lagrangian formalism, both sets of equations ( p.6D -( prTUD 
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and ( p.l2| )-( p.l4| ) are equivalent to the following second-order Langevin equation for 
the field (p: 

4>n + l<Pn = ^ + Vn (3.19) 



The second term in the l.h.s. of Eqn. ( p.l9| ) is the damping force. It dominates over 
the first one at large friction or large times t > 7"^ . Then neglecting the first term 
in the l.h.s. of Eq.( |3.19| ) one obtains the first order Langevin equation: 



M^ + 1) - M^) = - - ^ + J^ a (3.20) 

where 6 = XT and ^ is a random variable with Gauss's distribution of variance 1. ^ is 
the only parameter of the classical model (apart from the lattice spacing). To fix the 
physical temperature one needs to know A. Since A is a parameter of the semiclassical 
expansion according to Eqn.( ^^ , it is fixed only in the underlying quantum theory. 
The r.h.s of Eqn. (|3. 2 Op depends on the ratio e/7, not separately on e and 7. 
This implies a simple friction dependence of, say, the diffusion rate determined by 
means of the first-order Langevin equation: the diffusion rate must be inversely 
proportional to the friction coefficient r(T, 7) ~ I/7. Thus within the first-order 
Langevin treatment the absolute value of the diffusion rate is fixed in a simple way by 
the friction coefficient 7, which is a dimensional quantity. The second-order Langevin 
equation, valid for arbitrary friction, predicts nontrivial friction dependence, which 
has been evaluated explicitly by Kramers ^ in the case of one degree of freedom 
coupled to the heat bath. In the limit of large 7 that friction dependence of the 
escape rate, of course, reduces to the one we just derived from the first-order Langevin 
equation. The question we would like to clarify is how large the friction coefficient 
must be in the case of field theory to allow one to use the first-order Langevin equation 
instead of the second-order one. Our numerical data proves to be helpful for this end. 



10 



4 Numerical simulations 



We have solved the second-order Langevin equations (|3.15| - |3T7|) numerically for the 



system of size L = 50, a = 1. This volume is sufficient to accomodate a few kinks 
(which size is 1 in our units) at low temperatures. The diffusion parameter A(t) has 
been measured, defined following p as: 

(A(t)) = i ['"df {^{t' + t)-^{t')}' (4.1) 

to Jo 

The system was observed over 5 ■ 10^ to 2 ■ 10^ iterations (i.e. a time to ~ 10^). 
The longer runs were necessary to reduce statistical errors at the lower temperatures 
where the diffusion rate is smaller. The time step was typically e = 0.02 . The 
Brownian motion of ip has been unambiguously observed. Fig. 2 demonstrates the 
diffusion at temperature T = 6 for 3 different values of the lattice spacing. The rate 
is obviously insensitive to the lattice spacing and the relatively large value a = 1 may 
be used to obtain physical results. Large values of a are preferable to save computer 
time. 

Periodic as well as free boundary conditions were studied. In both cases we 
measured the diffusion rate in the range of temperatures 6 = 1.33 4- 18. The 
semiclassical domain corresponds to ^ 8 in accordance with ( |3.2| ). Fig. 3 shows 
the temperature dependence of the diffusion rate on a logarithmic scale. In the case 
of free boundary conditions all the data points starting from T = 3 downwards lie 
on a straight line. This implies the expected thermal activation behaviour (|3.5| ). 
The slope of the straight line is seen to be 8, which is exactly the classical value 
of the kink mass ( |3.2| ). In the case of periodic boundary conditions, a somewhat 
more interesting phenomenon is observed: starting around T = 2 the slope of the 
straight line changes from 8 to 16 , which corresponds to the classical energy of 
a kink-antikink pair. This is again in accordance with the expectations (|3.4| ) for 
periodic boundary conditions. The sensitivity to the boundary conditions appears in 
the lowest temperature domain. Fig. 4 shows that at higher temperatures there is no 
dependence on the boundary conditions. The cross-over temperature below which 
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the kink- ant ikink pair energy is observed in ( p.4| ) depends on the system size L. The 
larger the volume L the smaller that crossover temperature. 

The diffusion rate can also be measured in the high temperature domain ^ > 8 , 
where the semiclassical approximation is not valid. Our numerical results are shown 
on Fig. 5 for both types of boundary conditions. The diffusion rate does not depend 
on the choice of boundary conditions as it is not associated with the production of 
kink-antikink pairs any more. 

At high temperatures the Langevin equation ( ^.201 ) can be integrated explicitly, 
because the noise term, whose contribution is proportional to vT, dominates over 
the regular force. Then the analytical predition for the rate is : 

r = 4 (4.2) 



The numerical measurements are in perfect agreement with Eqn. (^4.2| ). 

The crossover temperature from the semiclassical to the high temperature domain 
is seen to be around T = 3. Perhaps surprisingly it is less than the kink energy (|3.2| ). 
To determine it more accurately we have measured the probability distribution of the 
averaged field (p, whose logarithm determines the effective potential V{ip^T). One 
can see on Fig. 6 that at T = 2 there is a deep parabolic well at the origin, so that 
the field most of the time oscillates around (^ = 0. Beyond the parabolic well the 
potential is fiat. This implies the free motion between the different vacuum sectors, 
which proceeds via the formation of kink-antikink pairs, seen as the modulation of 
the high frequency oscillations corresponding to the parabolic well. The deeper the 
well the better is the very notion of the kink configuration at T ^ 0. At temperature 
T = 3 the parabolic well is rather shallow, so the low frequency modulations usually 
interpreted as kinks become hard to identify. The effective potential is almost entirely 
fiat, which means that the transitions between different classical vacua do not require 
smooth low frequency configurations any more. This is obviously the boundary of 
the semiclassical domain of temperatures. At temperature T = 6 there is not any 
sign of the confining parabolic well. 
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Thus we have found that the effective potential V{(p, T) of the diffusing variable is 
flat in the high temperature domain, starting from T ~ 3. This result is interesting 
in the context of the baryogenesis within the standard electroweak theory [|l^ . In the 
Standard Model there is a variable, which is expected to diffuse at finite temperature. 
It is the Chern-Simons variable Nqs- Diffusion of Ncs leads to the dissipation of the 
baryonic number. This effect could be responsible for the production of the baryon 
asymmetry in the expanding Universe. In |jl^ a scenario for baryogenesis within the 



Standard Model was suggested, in which the crucial assumption was the fiat effec- 
tive potential V{Ncs) in fhe high temperature domain. While to measure properly 
V{Ncs) in fhe Electroweak Theory is not easy, our data in the Sine-Gordon model 
indicates that the fiat effective potential V{Ncs) is a rather plausible assumption. 
We would like to emphasize again the similarity between the Sine-Gordon model and 
Abelian Higgs theory in 1+1 - dimensions. The zero mode if of the Sine-Gordon 
theory plays the role of the Chern-Simons variable in the Abelian Higgs model. One 
can decompose </) ^ (^ + 5ip, to demonstrate that the effective potential V{ip,T) is 
a bounded periodic function, allowing for the Brownian motion. 

We now turn to the dependence on the friction coefficient. The friction coefficient 
7 is a dimensional input parameter in the Langevin equation ( p.l9|) . It determines 



the scale as it is seen, for instance, from formula ( [4.2| ), derived for large 7. Although 
the friction dependence of the escape rate has been obtained analytically in |]I| in 
the case of one degree of freedom exposed to the heat bath, it is a less trivial task 
in the field theory. We have measured the friction dependence of the diffusion rate 
numerically by solving equations (|3.15|) - ( p.l6|) for the two temperatures T = 2, 6 



and for various friction coefficients. The results are shown in Fig. 7, 8. In both 
the semiclassical and high temperature domains we find the same behaviour of the 
rate. The inverse rate varies linearly with the friction coefficient from large 7 > 1 
down to very small 7 ~ 10"^. This simple friction dependence coincides with the 
one predicted by the first order Langevin equation. Therefore we obtain experimental 
evidence that the first order Langevin equation may be used to explore the large time 

13 



behaviour of the correlators not only for large, but also for small friction coefficients, 
in particular, for 7 ~ 1. 

One can see from Fig. 8 that the rate does not diverge in the limit 7 ^ 0. It 
approaches some finite asymptotic value. Meanwhile the original Langevin equations 
become the conventional Hamiltonian ones in the limit 7-^0. Therefore we expect 
the rate extracted from the second order Langevin simulations in the limit of vanishing 
friction to coincide with the one obtained by means of microcanonical simulations. 
To verify this conjecture we first of all notice that the microcanonical simulations 
of classical systems in real time [1^ naturally correspond to some fixed energy, not 



temperature. To obtain the time dependent Gibbs averages, one has to average the 
result of the microcanonical measurement over the initial field configuration with the 
Gibbs weight: 

/A/,xv ^ / Vi^oVpo exp ( - ^(^0, Po) /T) {<^{t, (fp , pp) - <^{0) }^ 

^ ^'^ IVy,oVpoexp{-H{^o,Po)/T) ^ ' > 

The problem here is that the value of the diffusion parameter A at time t, ob- 
tained as a result of the Hamiltonian evolution, is a functional of the initial field 
configuration {fo,Po} = {fit = 0),p(t = 0)}. This functional is not known ex- 



plicitly. However, following the suggestion by J.Smit [|1^] we evaluate this functional 
numerically. This means that to measure the Gibbs averaged diffusion parameter 
directly we calculate numerically {f{t) — v^(0))^ from the Hamiltonian evolution for 
a given starting configuration {(po,Po}, then update this starting configuration by 
means of the Metropolis procedure, corresponding to some temperature T, then do 
microcanonical evolution again and so perform the Gibbs average. 

In this way the diffusion parameter (|4.3| ) has been measured. The function (A(t)) 
is shown on Fig. 9. One can see a crossover between two domains of time. At short 
times t < O(IO^) the deterministic behaviour is observed (A(t)) ~ t^, while at 
t > O(IO^) the Brownian motion sets in (A(t)) ~ t. The corresponding diffusion 
rate as a function of temperature is shown on Fig. 10. The thermal activation 
behaviour is clearly seen at low temperatures. As a result we observe a remarkable 
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coincidence between the data from the second order Langevin simulations in the hmit 
of vanishing friction and the direct Gibbs averages. One obtains r( T = 2 ) ~ .3 and 
r(T = 6) ^ 8. from both Fig. 8 and Fig. 10. This resuh answers the question of 
how to relate the Langevin measurements to the direct calculations of the Gibbs's 
average ( |4.3D . 



5 Comments on the errors 

As mentioned above all the runs which results are presented were long enough to 
make statistical errors negligible. The two main sources of systematic errors are the 
finite lattice spacing a and the time step e. The insensitivity to the lattice spacing 
was discussed before (see Fig. 2). Here we would like to address the artifacts of the 
discretization of the Langevin time. 



It is known |T9| that the finite time step in Langevin simulations makes the tem- 
perature of the simulated system Tg// larger than the input one T: AT = Teff — T ~ 
0{e). To check the effective temperature we measure the diffusion parameter at rather 
small times, as shown on Fig. 11. The Brownian diffusion is clearly observed at short 
times, but it has nothing to do with the motion of the system between topologically 
different vacua. At short times the white noise always dominates over the regular 
force in the Langevin equation. Therefore the short time behaviour of (A(t)) is, in 
fact, controlled by the free Langevin equation, which immediately implies diffusion 
of 1^. This short time diffusion is already sensitive to the discretization effects. Its 
diffusion rate is given by Tg/j, which has been measured numerically and is shown on 
Fig. 12. One can see that the effective temperature of the simulated system follows 
the input one with fairly high accuracy. Since these measurements are obtained from 
very short runs, the computer time needed is quite small. So we have also checked 
that AT = ce, with c = 1.5 ± .5. 

Another alternative to check the effective temperature of the simulated system is 
to measure the Langevin average of the canonical momentum squared over the whole 
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run: < p^ > = T^ff. This way is not easily generalizable for theories where there is 
a couphng between the canonical momenta and coordinates in the Hamiltonian, like 
in gauge theories. 

6 Conclusions 

The study of the Sine-Gordon field theory in 1+1-dimensions shows that the classical 
Langevin simulations in real time prove to be efficient in obtaining valuable informa- 
tion about the nonperturbative effects in field theory at high temperatures. In the 
semiclassical domain of temperatures kinks are seen as smooth modulations of the 
high frequency oscillations of the field. Diffusion between different classical vacua is 
due to random process of production and free motion of kink-antikink pairs during 
the Langevin evolution. This is very well confirmed by the measured temperature 
dependence of the diffusion rate, which exhibits thermal activation behaviour with 
the classical value of the kink mass in the Boltzmann exponent. 

Crossover is observed between the semiclassical and high-temperature domain, 
where the semiclassical approximation is not valid. The effective potential of the 
diffusing variable is found fiat at high temperatures. In the high temperature domain 
at moderate friction the diffusion rate follows prediction of the free Langevin equation. 

The dependence of the diffusion rate on the friction coefficient is found identical 
to the prediction of the first order Langevin equation down to rather small friction 
coefficients 7 ~ 10^^. This justifies the use of the first order Langevin equation in 
the calculations of correlation functions at large real times. 

We also measure the diffusion parameter directly as the Gibbs ensemble average by 
means of the Metropolis procedure and microcanonical evolution. The correspond- 
ing diffusion rate is found to coincide with the smooth extrapolation of Langevin 
measurements in the limit of vanishing friction coefficient. 
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Figure Captions 

Fig. 1 Typical potential creating a metastable state at finite temperature. 

Fig. 2 Time dependence of the diffusion parameter (O!) obtained from the second 



order Langevin equation. Measurements were performed every 100 iterations. The 
parameters are: T = 6,L = 50, e = .02. Three solid lines correspond to different 
lattice spacings: a = {1, .5, .25} ; the dashed line is a fit with slope 1. 

Fig. 3 Temperature dependence of the diffusion rate on a logarithmic scale for: 
a) free, b) periodic boundary conditions. L = 50, a = 1. Typical uncertainty is 
shown for T = 2. The dashed lines correspond to the Boltzmann exponent with 
slope given by the classical value of the kink mass M^ = 8 in a) and the energy of 
the kink-antikink pair 2Mfc = 16 in b). 

Fig. 4 Temperature dependence for free vs. periodic boundary conditions at both 
high and low temperatures. 

Fig. 5 Temperature dependence of the diffusion rate at high temperatures. The 
dashed line is a prediction from the corresponding first order free Langevin equation. 

Fig. 6 Logarithm of the probability distribution (effective potential) of (f at 
different temperatures indicated on the plots. 

Fig. 7 Friction dependence of the rate in the low (a) and high (6) temperature 
domains obtained from the second order Langevin equation. L = 50, a = 1. The 
dashed straight line is a prediction from the first order Langevin equation, derived 
normally for large friction 7 > 1 or large times. 

Fig. 8 Same as in Fig. 7, but only the points corresponding to very small values 
of the friction coefficient are shown. The cross-points on the y-axis are the extrapo- 
lations. 

Fig. 9 Time dependence of the diffusion parameter ( [4.1D obtained from direct mea- 



surements of the Gibbs average (43) for periodic boundary conditions. Temperature 
is T = 1.5. The initial field configuration was updated by means of the Metropo- 
lis procedure 100 times and each time 10^ leap-frog iterations of the microcanonical 
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evolution were performed. L = 50, a = 1, e = .1 

Fig. 10 Temperature dependence of the diffusion rate obtained from direct mea- 
surements of the Gibbs average ( [4.3|) . The dashed hne corresponds to the Boltzmann 
factor with the energy of the kink-antikink pair. The parameters of the measurements 
are as in Fig. 9. 

Fig. 11 The typical diffusion at short times, controlled by the free Langevin 
equation. 

Fig. 12 The effective temperature of the simulated system obtained from the 
short time diffusion. The dashed line is a fit. The solid line corresponds to the input 
temperature. The deviation is proportional to the time step e. 
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